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Our experimental studies of electron transport in wide (14 nm) HgTe quantum wells confirm 
persistence of a two-dimensional topological insulator state reported previously for narrower wells, 
where it was justified theoretically. Comparison of local and nonlocal resistance measurements indi¬ 
cate edge state transport in the samples of about 1 mm size at temperatures below 1 K. Temperature 
dependence of the resistances suggests an insulating gap of the order of a few meV. In samples with 
sizes smaller than 10 /rm a quasiballistic transport via the edge states is observed. 


The topological insulators (TI) represent a quantum 
state of condensed matter with insulating bulk and con¬ 
ducting gapless states at the surface or edge [H-Q- The 
existence of such materials is justified within a concept of 
topological ordering introducing order parameters which 
are often expressed as invariant integrals over the mo¬ 
mentum space. In the presence of time reversal symme¬ 
try, the materials with energy band gaps (band insula¬ 
tors) are classified by Z 2 topological invariants @ which 
take two values, 1 or 0, thereby providing a distinction 
between topological and normal insulators. Mathemati¬ 
cally, one can construct Z 2 invariants in different ways, 
but their physical meaning always relies on the symmetry 
of electron wave function which is changed as a result of 
energy band inversion. Such an inversion occurs due to 
spin-orbit coupling and Darwin term contributions in the 
Hamiltonians of the crystals formed from heavy atoms. 
There are three types of band inversions {s — p, p — p, 
and d — f) in the three-dimensional (3D) TI discovered 
so far 

The most extensively studied TI materials, bismuth 
chalcogenides and related alloys, belong to p—p inversion 
type. Eor thin layers of these materials, one expects a di¬ 
mensionality crossover: when layer thickness d decreases, 
the material transforms from 3D TI into two-dimensional 
(2D) TI 0. This occurs when the wave function de¬ 
cay length of the surface states becomes comparable to 
d. As a result, the 2D states from opposite surfaces hy¬ 
bridize and their spectrum is no longer gapless 0, @ ■ On 
the other hand, since the surface states in TI cover the 
whole surface of the layer, including side regions, they are 
transformed into ID conducting edge channels in these 
regions, Fig. 1 (a), 

A special place in this connection belongs to HgTe, a 
zinc-blend-structure crystal with s — p band inversion, 
where the energy of p-type Tg band in the T point of 
Brillouin zone is higher than the energy of s-type Tg 
band. In spite of the inversion, bulk HgTe is not a 3D 
TI because it is a symmetry-protected gapless material. 
A gap can be opened in a thin HgTe layer sandwiched 
between Cdi_a;Hg 3 ;Te layers (normal insulators), as re¬ 
alized in epitaxially grown quantum wells (QWs). Due 


to size quantization, the heavy hole (hh) continuum that 
forms the valence band in HgTe splits into a set of 2D 
states hybridized with two interface-like states S'! and 
S2, Fig. 1 (b). The gap between the ground-state hh 
subband (hhl) and the next adjacent subband exists for 
QW narrower than 18 nm (for wider wells the QW is in 
a semimetallic state). Since the gap opening is accom¬ 
panied with the dimensionality crossover, a HgTe-based 
QW should be a 2D TI having edge states in the gap 
between subbands. However, a direct theoretical proof 
for this statement, based on a two-subband effective 2D 
Hamiltonian 0 , has been done only for a special situa¬ 
tion when 51 subband is just slightly below hhl subband, 
which is applicable to narrow wells in the width range of 
approximately 6.3 — 8.3 nm. The edge state transport 
in such QWs was also confirmed experimentally m- In 
wider HgTe QWs, any effective Hamiltonian methods are 
not generally feasible because of a complicated subband 
structure, though a usage of three-subband (51, hhl, and 
hh2) basis [l^ extends the range of applicability of such 
methods. In particular, it was found that when 51 sub¬ 
band falls below hh2 one (so the principal gap is formed 
between hhl and hh2 subbands) the edge states exist 
both in this gap and in the next gap between hh2 and 
51 subbands. 

A question that naturally arises concerns persistence of 
2D TI state in wide QWs where the 51 subband lies be¬ 
low several hh subbands, so the situation is far different 
from that described theoretically in Refs. [13 and [T3 - 
From the point of view of dimensionality crossover, there 
are no reasons to deny the 2D TI nature of these sys¬ 
tems, since widening of the QW (actually, approaching 
of HgTe layer to the bulk state) does not cancel the fact 
of s — p inversion in this layer and, accordingly, cannot 
destroy the edge states. In this Letter we report exper¬ 
imental investigation of 14 nm wide HgTe QWs which 
are wider than those studied previously [ll|, O, [g but 
still have a sizeable gap of a few meV between hhl and 
hh2 subbands. We plot the expected edge states in this 
gap schematically in Fig. 1 (c) as two (one for each spin 
number) gapless branches merging with 2D subbands on 
a tangent. By transport measurements, we indeed ob- 
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FIG. 1: (Color online) (a) Cross-section of a 3D TI sample 
(schematic) with surface state channels shown in red. As the 
thickness of the sample decreases, a transition to 2D TI takes 
place, (b) Energy spectrum of size-quantized subbands in a 
symmetric 14 nm wide Cdo. 65 Hgo. 35 Te/HgTe/Cdo. 65 Hgo. 35 Te 
QW, calculated numerically by using Kane Hamiltonian. 
Dashed (blue) lines show the spectrum of interface states at a 
single HgTe/Cdo. 65 Hgo, 35 Te boundary under approximation 
that mixing of these states with hh states is neglected, (c) A 
magnified picture of two upper hh subbands and expected 
edge state spectrum, thin (red) lines, in the gap between 
them. 


tain numerous proofs for the edge state transport in these 
QWs. 

The experimental samples are Hall bridges fabricated 
on top of the 14 nm wide HgTe QW with the surface 
orientation (112) and provided with electrostatic gate . 
Their fabrication technology is described in detail in . 
Three different types of experimental samples were used: 
macroscopic Hall bridges (see Fig. 2) with the width 
50 /rm and the distance between the voltage probes 100 
fim and 250 /im, and two types of microscopic samples, 
whose layout together with the scale is shown in Fig. 3 
M- The transport measurements were conducted in the 
temperature range 0.2 — 10 K and in magnetic fields up 
to 10 T using the standard phase detecting scheme on 



E 1000 
6 

500 


(b) A 

T=0.2 K 

T=0.3 K 


T=0.4 K 

J 1 

T=0.6 K 

R /A\ 1 

T=0.8 K 

l=4-8;V=5-9 / / \ 1 

T=1 K 

/ / \ 1 

T=1.5K 

_ 

T=3.3 K 


-4 


V (V) 


-2 


E 

■W ,»» » * 

0 

»» ® ® ® 

10^ 


® * * 


• « 1 10- 



~ s’” 

■a Calculations 


® R, ^,0. 

' a R, 

10 ' 

’® ® R 

® « R^, 


nl 

2 ifr(K'') * 


1/T(K'’) 


FIG. 2: (Color online) Local (a) and nonlocal (b) resistance 
at different temperatures, sample I, B = 0 (insert - schematic 
view of the sample). (c) Temperature dependence of local and 
nonlocal resistance when the Fermi energy is situated in the 
center of the insulating bulk gap (insert - calculation). 


frequencies 3 — 12 Hz and the driving current 0.1 — 1 nA 
to avoid heating effects. The electron mobility /r in all 
samples studied was above 10® cm^/V s for the carrier 
density 3 x 10^^ cm“^. 

We first consider the properties of our macroscopic 
samples. Observation of both local (i?/) and nonlocal 
(Rni) resistance in the absence of magnetic field is gener¬ 
ally considered as a direct proof of 2D TI state. In that 
case the current flows through a sample along its borders 
and resistance exists regardless to positions of voltage 
probes with respect to current contacts, either in ballistic 
[T^ or diffusive [3 regime. Figures 2 (a) and 2 (b) show 
Ri and as functions of the gate voltage for various 
temperatures in the range 0.2 —3.3 K and Fig. 2 (c) shows 
these resistances as functions of temperature for the gate 
voltages corresponding to the maxima of local and non¬ 
local resistance, i.e. when the Fermi level is in the middle 
of the insulating gap. Qualitatively, the behavior of local 
and nonlocal resistances is similar. When the temper¬ 
ature is above 1 K, both resistances grow exponentially 
with decreasing temperature. For lower temperatures the 
resistance growth slows down and for T < 0.5 K it obeys 
a power law R oc T”" {a ~ 0.5). As this takes place, 
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even though at higher temperatures the local resistance 
is several orders of magnitude greater than the nonlocal 
one, at T < 1 K the nonlocal resistance becomes larger. 
The described behavior of Ri and Rni is typical for the 
2D TI, as it follows from the fundamental difference in 
the relative contributions of the bulk and edge trans¬ 
port when measuring in local and nonlocal configurations 
[l4l| . At higher temperatures, when the bulk contribution 
is still sufficiently large, Rni is exponentially small com¬ 
pared to the sheet resistivity. With decreasing T the bulk 
contribution to transport also decreases, and at some T, 
depending on the insulating bulk gap, becomes negligi¬ 
bly small. Under these conditions the difference in the 
resistance values measured in local and nonlocal configu¬ 
rations is determined only by the distribution of currents 
flowing along the sample perimeter and by the position of 
the voltage probes. At T < 0.5 K for all investigated con¬ 
figurations the resistance at its maximum is more than 
an order of magnitude greater than /i/2e^ which means 
that the transport via the edge states is diffusive. 

Let us discuss in more detail the temperature depen¬ 
dence of Rni and Ri. First, the maximum of the curves 
Rni{Vg) and Ri{Vg) shifts to the right with the temper¬ 
ature decreasing. Such behavior was not observed in 
2D TI previously and is probably related to the com¬ 
plicated energy spectrum of the system investigated, or, 
more specifically to its much smaller gap which is fur¬ 
ther diminished by the bulk bands density of states tails. 
Second, as has already been mentioned, with lowering 
temperature the resistance, after the exponential growth, 
continues to increase, but at a much lower rate, following 
a power law R oc typical for quasi ID wires in a 

weak localization regime. We have attempted a quanti¬ 
tative description of the temperature dependence of the 
local and nonlocal resistance peak values by using the 
model proposed in ©[II. The activation energy for 
the bulk transport has been chosen as a fitting param¬ 
eter. The results of the calculation presented in Fig. 2 
(c) show a reasonable agreement between the calculated 
and measured dependences. The value of the activation 
energy found from fitting the calculation to experiment 
is approximately 1.2 meV for local and nonlocal configu¬ 
rations alike, as was expected. This value is less than A= 
3.3 meV, the indirect insulating gap value obtained from 
the energy spectrum calculation. This discrepancy is not 
surprising if one considers the disorder due to impurities 
and QW thickness fluctuations that is always present in 
a real HgTe sample. In such case one may expect that 
the activation energy would correspond to the mobility 
gap rather than to a much larger calculated insulating 
bulk gap ©I . 

The mean free path for the edge states transport de¬ 
termined experimentally for the diffusive transport in the 
macroscopic samples is 2 — 5 /rm for sample I and 12 — 14 
/im for sample II (See supplementary material). On this 
account our QWs look promising for observation of bal¬ 
listic transport via edge states, considering that in the 
majority of the previously studied 2D TI this value was 


close to 1 fim. For this purpose, two types of micro¬ 
scopic samples were fabricated, one with the dimensions 
W = 1.7 ^m and L — 1.8 ^m (see Insert to Fig.3a) and 
the other an H-shaped bridge with the width W = 3.2 
fim and the length L = 2.8 p,m (see Insert to Fig.3b). 
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FIG. 3: (Color online) R{Vg) dependences for local (red 
curves) and nonlocal (green curves) resistance measurement 
configurations obtained in two different types of microscopic 
samples. The samples layouts are shown schematically in the 
Inserts. All curves were obtained at 0.2 K except the black 
curves measured at temperatures 4.2 and 4.6 K. Opposite 
each curve the corresponding measurement configuration is 
shown schematically. The dashed horizontal lines mark the 
resistance values expected for these configurations in the case 
of a purely ballistic transport via edge current states. 

Figure 3 shows the dependence of the local (red color) 
and nonlocal (green color) resistance versus gate voltage 
measured in the two microscopic samples. When the bias 
applied to the gate decreases below —IV (Fig.3a) (2 V 
in Fig.3b), the local resistance starts to grow gradually 
as the Fermi level first descends to the bottom of the 
conduction band and then enters the insulating bulk gap 
at Vg w -4.5 V (Fig.3a) {Vg « -2 V in Fig.3b). The 
nonlocal response is close to zero when the Fermi level 
remains in the conduction band. The abrupt increase in 
the nonlocal resistance as the Fermi level enters the insu¬ 
lating gap signals, as well as in the macroscopic case, the 
presence of the edge current states in the gap. The calcu¬ 
lation of the local and nonlocal resistance values expected 
in our microstructures in the case of ballistic transport 
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via edge states is quite simple and is indicated by the 
dashed lines for each configuration shown in Fig.3. As 
may be seen, the average resistance values measured at 
T = 0.2 K in the gate voltage range corresponding to 
the Fermi level staying in the insulating gap are quite 
close to the levels expected for purely ballistic transport. 
With the negative bias increasing the Fermi level enters 
the valence band and the sample resistance decreases. 
Relying on this data, it is possible to conclude that our 
microstructures demonstrate a quasiballistic edge trans¬ 
port in a 2D TI. This fact is quite important considering 
that up till now the observations of ballistic transport in 
HgTe-based 2D TI reported in m has remained unique. 
We also observe other similarities with the results re¬ 
ported in [HI- In particular, when the Fermi level lies 
in the insulating gap there are random fluctuations both 
in the local and nonlocal resistance. The amplitude of 
these fluctuations sharply decreases as the temperature 
increases indicating their mesoscopic origin. 

Figure 3 also shows the variation of the local resistance 
with temperature in the range from 0.2 K to 4.2K in our 
microstructures. One can see that this variation is notice¬ 
ably weaker than in our macroscopic samples (compare 
with Fig. 2 (c)). This fact has a simple explanation. In 
microstructures the sheet conductance due to the bulk 
transport is the same as in macroscopic samples while 
the resistance to transport via edge states decreases by 
approximately one order of magnitude. This observa¬ 


tion also indicates that the mean free path for the edge 
states transport in our QWs must be comparable to or 
higher than « 10 /rm. We also observe a suppression 
of both local and nonlocal quasiballistic conductance by 
weak magnetic field, a behavior typical for 2D TI state 
[ullii , which will be reported elsewhere. 

We believe that our observation of relatively high mean 
free path for the edge states is not accidental but rather 
related to the advantages associated with the use of a 
wider quantum well. Indeed, the width of any QW is not 
uniform but fluctuates from point to point around its av¬ 
erage value d with the amplitude S. That S is determined 
by the growth technology employed and is practically in¬ 
dependent of the QW width. The fluctuation of the QW 
width results in a random potential in the bulk of the 
QW. However, the amplitude of that random potential 
would be much smaller in a wider QW as it is propor¬ 
tional to l/{d^). From that point of view it is clear that a 
wider QW well is more advantageous for the observation 
of ballistic transport in 2D TI. 

The results of the present study confirm that the 2D TI 
state in HgTe QWs is quite robust and exists in a sizeable 
range of well widths despite of the fact that the energy 
spectrum in such QWs is complicated and strongly de¬ 
pendent on the well width. 

The work was supported by the RFBI grants, by the 
RAS grant 24.11 and by FAPESP, CNPq (Brazilian agen¬ 
cies). 
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In this supplementary, we provide details on the sample description, the mean free path calculation 
and the model taking into account the edge and bulk contribution to the total current. 


SAMPLE DESCRIPTION 

In the main text we describe the experimental results 
in macroscopic and mesoscopic samples in zero magnetic 
field. The Cdo.esHgo.asTe/HgTe/Cdo.esHgo.asTe quan¬ 
tum wells with surface orientation (II2) studied here 
were grown by molecular beam epitaxy (MBE). The 
schematic view of the samples layer structure is shown 
in Fig. la. As shown in previous publications [l|, the 
use of substrates inclined to the singular orientations 
may lead to the growth of more perfect films. There¬ 
fore, the efficient growth of alloys is done predominantly 
on substrates with surface orientations of (112), which 
deviate from the singular planes by approximately 19°. 
We have investigated samples of both macroscopic and 
microscopic dimensions. The microscopic devices were of 
two types. The first with the dimensions W = 1.73 ^m 
and L = 1.8 /rm and the other an H-type bridge with 
the width W = 3.2 /rm and the length L = 2.8 ^m. Fig. 
1 (b). The macroscopic sample consists of 3 narrow (50 
/rm wide) consecutive segments of different length (100, 
250, 100 irm), and 8 voltage probes. Fig. 1(c). The ohmic 
contacts to the two-dimensional gas were formed by the 
in-burning of indium. To prepare the gate, a dielectric 
layer containing 100 nm Si 02 and 100 nm Si 3 Ni 4 was 
first grown on the structure using the plasmochemical 
method. Then, the TiAu gate was deposited. 

Low-temperature magnetotransport measurements 
were carried out to identify the charge neutrality point 
(CNP) and to obtain the relation between the induced 
carrier density and the gate voltage Vg. Sweeping the 
gate voltage, we have found a continuous crossover from 
electron to hole type of conductivity, when the Fermi 
level is shifted from the conduction to the valence band. 
The possibility to tune the carrier density electrostati¬ 
cally enables the observation of electron and hole trans¬ 
port in a single device, making it possible to control the 
filling of Landau levels in the presence of the perpendic¬ 
ular magnetic field. Figure 2 shows longitudinal p^x and 
Hall pxy resistivities as functions of the gate voltage at 
a fixed magnetic filed H = 2 T. Pronounced plateaux 
with values pxy = h/ve^ are clearly seen for electrons at 




FIG. 1: (Color online) (a) Schematic layer strncture of a wafer 
containing 14 nm HgTe QW. (b) mesoscopic samples topol¬ 
ogy, yellow rectangnlar area depicts the gate; (c) macroscopic 
samples layout. 


V from —1 to —5, accompanied by deep minima in Rxx- 
As Vg is swept through the CNP, when the Fermi level 
is in the middle of the gap, the longitudinal resistivity 
shows a maximum, whereas pxy goes gradually through 
zero from negative on the electron side to positive on the 
hole side. It is a clear indication of CNP position and 
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electron-hole crossover as a function of the gate voltage. 
Observation of the edge state transport near CNP in zero 
magnetic field is described in the main text. From the 
resistance minima we obtain the density variation with 
gate voltage 8.4 x 10^'* which coincides with the 

density extracted from the capacitance. 
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TABLE I: The first column: sample number. The second 
column: various local and nonlocal measurement configura¬ 
tions. The third column: mean free path for the edge states 
transport obtained from the comparison of the peak resistance 
value with theory. 

On the other hand, it is in agreement with a conclusion 
made earlier ii concerning the absence of topologi¬ 
cal protection against backscattering on the scale of few 
microns and longer. 


FIG. 2: (Color online) (a) Longitudinal (black line) and 
Hall pxy (red line) resistivities as functions of the gate voltage 
at a fixed magnetic field, B = 2.0 T. Vertical arrow indicates 
CNP. 


MEAN FREE PATH CALCULATION FOR THE 
EDGE STATES TRANSPORT 

At T < 0.5 K for all investigated measurement con¬ 
figurations in the macroscopic sample the resistance at 
its maximum is more than an order of magnitude greater 
than h/2e^ which means that the transport via the edge 
states is diffusive. In this case there is a simple ex¬ 
pression that allows one to calculate the resistance value 
for any measurement configuration assuming that there 
is only edge state transport in the sample: = 

where is the resistance measured 

between contacts i and j while the current is maintained 
between contacts n and m, Li j {Ln^m) are the distances 
between i and j {n and m) along the sample edge that 
does not include n and m (i and j), L is the total perime¬ 
ter of the sample, and I is the mean free path. Thus, it is 
possible to determine the mean free path I for the trans¬ 
port via the edge states. Table 1 presents the results 
obtained in this way for samples I and II. The values of I 
derived from different measurement configurations in the 
same sample are quite close, but may differ substantially 
in different samples. This data supports the statement 
advanced in Q that in 2D TI fabricated on the basis 
of HgTe QWs with inverted energy spectrum the edge 
states have a macroscopic length (about a millimeter). 


EDGE+BULK MODEL 

In the main text it has been demonstrated that the 
resistance of HgTe quantum wells reveals a broad peak 
when the gate voltage induces an additional charge den¬ 
sity altering the quantum wells from n-type conductor 
to p-type one around the charge neutrality point (CNP). 
The peak amplitude is approximately equal to /i/2e^ in 
mesoscopic samples, due to nearly ballistic edge-state 
transport, and R S> /i/2e^ in macroscopic samples. The 
lack of the robustness of resistance quantization against 
intrinsic and introduced disorder in macroscopic samples 
may be explained in various ways Independently 

of the particular microscopic mechanism responsible for 
this, the resistance can be described by the combination 
of the edge-state and bulk transport contributions, with 
taking into account both the backscattering within one 
edge and bulk-edge coupling. The local and nonlocal 
transport coefficients arise from the edge state contri¬ 
bution and short-circuiting of the edge transport by the 
bulk contribution, the latter is more important away from 
the CNP. Below we reproduce the basic features of this 
model [ 3,0 . The transport properties in the bulk can be 
described by the current-potential relation 

tW = (1) 


where i = 1,2 numbers the states with different projec¬ 
tions of the spin, ijji are the electrochemical potential for 
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electrons, and r = {x, y) is the 2D coordinate. As we 
consider the isotropic conduction, the non-diagonal part 
of the conductivity tensor appears only at non-zero mag¬ 
netic field. Assuming the components of the conductiv¬ 
ity tensor as coordinate-independent parameters, we can 
solve the problem by solving the Laplace equation for the 
potentials, = 0 , because the charge conservation 

continuity conditions requires Vjj(r) = 0. A solution to 
Laplace equation is uniquely determined by specifying 
boundary conditions, which in our case are modified by 
the bulk-edge current leakage. In order to describe the 
transport in the presence of the edge states, we introduce 
two phenomenological constants 7 and g, which represent 
edge to edge and bulk to edge inverse scattering length, 
respectively. Then, the boundary conditions expressing 
zero current normal to the boundary in the presence of 
the bulk-edge coupling are given by 

nj, = gi^Pi - (2) 


For both local and non-local configurations, the resis¬ 
tance is calculated as 

Rxx — ^tot 5 ^tot — ^edge “t” ^bulk: 


^ = 2 (‘^11 “+ *^21 - V521'). (5) 


where V is the potential difference at the voltage probes, 
Itot is the total current flowing between current contacts, 
ifii and ifiy are the potentials at the voltage probe loca¬ 
tions. The edge and bulk currents for the local case (at 
arbitrary point x along the sample) are given by 


ledge = “ <^2 + ^2' - : 


ibulk 


= E 


^xv - O^xl 


■ diiii 
OX 


( 6 ) 


where ipi are the local chemical potentials of the edge 
states, ijji 3'iid are the potentials and currents at the 
boundary, and n is a unit vector normal to the boundary. 

The edge state transport can be described by conti¬ 
nuity equations taking into account the scattering 

between edge and bulk: 

= 7(¥’2 - +5(V'i - V^i), (3) 

-dxg}2=i{vi -‘P 2 ) + g{.'ip2-V2)- (4) 

The general solution of this problem, therefore, in¬ 
cludes solution of a 2D Laplace equation for the bulk 
electrochemical potentials tjji, 2 {x,y) together with Eqs. 
(2),(3),(4) describing the scattering between edge states 
and between edge and bulk states. The current can be 
calculated from this solution as a sum of the contribu¬ 
tions from bulk and edge states. 

In nonlocal configurations the edge -|- bulk model can 
be solved only numerically. We have performed self- 
consistent calculations to find Tpi ^2 solution of Laplace 
equation in two space dimensions and ^ 31^2 solutions of 
equations 3 and 4 on the edge using the Hall bar ge¬ 
ometry schematically shown in Fig. 1 (c) [or Fig. 2 of 
the main text]. The contacts are assumed to be thermal 
reservoirs, where full mixing of electron spin states and 
bulk states occurs Note that, in contrast to the stan¬ 
dard QHE, when mixing of the edge states occurs within 
metallic Ohmic contacts, in our samples potential mix¬ 
ing is provided by 2D electron gas in the region outside 
of the metallic gate. 

The equations for are discretized by the finite ele¬ 
ment method. The generalized Neumann boundary con¬ 
ditions, Eq. (2), are set in the regions outside the metal 
contacts. To solve the boundary value problem for a sys¬ 
tem of ordinary differential equations (3) and (4) we use 
a finite difference code that implements the 3-stage Lo- 
batto Ilia formula. The boundary conditions inside the 
metal contacts are set to ipi ^2 = ^ 1 , 2 - 


where (pi, 'tpi and (pii, ipii are the potentials at the oppo¬ 
site (bottom and top, respectively) edges of the sample, 
and the integral is taken across the sample from bottom 
to top. For non-local case the currents are calculated 
from similar expressions: 


ledge = “ ^^2 + ^2' - ^v) ■ 


^bulk — 


2 = 1,2 ^ 


o-yi (V'i - V'iO 


- cr^' 

^ XX 


[ dx^ 

dy. 


(7) 


where now (pi, ipi and pi', ipi' are the potentials at the 
opposite (left and right, respectively) edges of the cur¬ 
rent contact, and the integral is taken across this current 
contact from left to right. 

The conductivities are calculated as ai]} = = 

e^prin -f yLpp)/2, where /i„ and n {pp and p) are elec¬ 
tron (hole) mobilities and densities. To find n and p, the 
bulk densities of states for electrons and holes are repre¬ 
sented by steps Dn^p = 47 rTO„_p//i^ with mn,p being the 
effective masses of electrons and holes. The energy gap 
separating electron and hole bands is Eg = 3.3 meV. In 
addition, the sharp band edges are smoothed according to 
a Gaussian law with broadening energies F^^p = ?i/r„_p, 
where r„_p = Pn,piTin,p/e. The following parameters have 
been used: = 80000 cm^/V s, Pp = 5000 cm^/V s, 

rUn = 0.024 niQ, rup — 0.15 toq, where toq is the free 
electron mass. The extended states are separated by a 
mobility gap 1 meV. Our calculations reproduce the 
temperature dependence of the local and nonlocal resis¬ 
tances, shown in Fig. 2 of the main text. The best 
agreement between the experiment and theory is reached 
for parameters 7 “^ = g~^ = 5 pm, which roughly agrees 
with the mean free path I represented in the Table I of 
the main text. Parameters 7 and g are assumed to be 
independent of temperature. 
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